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Abstract 

Ph Bismuth ferrite, BiFeO^, is the only known room-temperature 'multiferroic' material. We 

o 

demonstrate here, using neutron scattering measurements in high quality single crystals, that the 
antiferromagnetic and ferroelectric orders are intimately coupled. Initially in a single ferroelectric 

"o 

c/3 state, our crystals have a canted antiferromagnetic structure describing a unique cycloid. Under 

electrical poling, polarisation re-orientation induces a spin flop. We argue here that the coupling 
between the two orders may be stronger in the bulk than that observed in thin films where the 
cycloid is absent. 
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Electricity and magnetism are properties which are closely linked to each other. This link 
is dynamic in essence, as moving charges generate a magnetic field and a changing magnetic 
field produces an electric field. This forms the basis of Maxwell's equations. In a solid, 
a similar coupling was first considered by Pierre Curie [1] between the magnetization M 
and electric polarization P. This magneto-electric (ME) effect was recently understood to 
be potentially important for applications because in information technology, it would allow 
magnetic information to be written electrically (with low energy consumption) and to be 
read magnetically. The ME effect was demonstrated and studied in the 1960s in Russia [2] 
and since then, many so called 'multiferroic' materials have been identified [3]. However, so 
far the magnitude and operating temperatures of any observed ME coupling have been too 
small for applications. In fact, the only known multiferroic material of potential practical 
interest is bismuth ferrite, BiFeO^ which is actually antiferromagnetic below T/v ~ 370 °C 
[2] and ferroelectric with a high Curie temperature: Tc 820 °C [4]. As a result, in recent 
years, there has been a resurgence in the research conducted on this material. Moreover, 
epitaxial strain in BiFe0 3 thin film has been described as a unique means of enhancing 
magnetic and ferroelectric properties [5]. It is actually unclear whether this is indeed the 
case and in order to clarify this point, the intrinsic properties of the bulk material need to 
be better understood. It is stunning that although BiFe0 3 has been extensively studied 
over the past 50 years, some of its most basic properties are still not fully known. For 
instance, it is only in 2007 that its spontaneous polarisation at room-temperature has been 
measured to be in excess of 100/iC/cm 2 [6]. Moreover, in the bulk, the coupling between 
magnetic and ferroelectric orders has never been fully clarified. This property has only been 
measured in thin films [7] very recently. This lack of accurate data stems from the difficulty 
in making high quality single crystals. We have recently been able to grow such single 
crystals below their ferroelectric Curie temperature using the flux technique [6, 8]. They 
are usually produced in the form of platelets 40-50 microns thick and up to 3mm 2 in area. 
Polarised light imaging and P(E) measurements [8] indicate that the as-grown crystals are 
generally in a single ferroelectric/ferroelastic domain state. We report here on a neutron 
study of the coupling between magnetic and ferroelectric orders in two of these crystals. 

Our BiFeOz single crystals are rhombohedral at room temperature with the space group 
R3c and a pseudo-cubic cell with a pc = 3.9581 A, a pc = 89.375°(a/ iea; = 5.567(8) A, 
Chex = 13.86(5) Ain the hexagonal setting), in perfect agreement with previous reported 
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data [9] . No ferroelastic twinning was observed and the elongated rhombohedral direction, 
which is parallel to the polarisation, is indexed as (111). Fe 3+ ions are ordered antiferro- 
magnetically (G-type) and their moments describe a cycloid with a period of 62nm, as has 
been established by early neutron diffraction data on sintered samples [10, 11]. Because of 
the rhombohedral symmetry, there are three equivalent propagation vectors for the cycloidal 
rotation: k\ = (50 — 5), k 2 = (0 — 55) and k 3 = (—550) where 5 = 0.0065. In powder neu- 

— * 

tron diffraction, the different equally populated k domains lead to a splitting of magnetic 
peaks along three directions. Thus, as has been pointed out recently [11], the determination 
of modulated magnetic ordering is not unique because elliptical cycloids and Spin Density 
Waves (SDW), give the same diffraction pattern. The exact nature of the periodic struc- 
ture is an important parameter for antiferromagnetic ferroelectrics since recent models of 
magnetoelectric coupling give a non-vanishing electric polarization for cycloids and elliptic 
ordering and zero polarization for a SDW [11, 12]. It is possible to eliminate this ambigu- 
ity by measuring high-resolution scans around the strongest magnetic reflections of a single 
crystal. This is however a difficult experimental challenge because the long period imposes 
an extremely high angular resolution. The diffractometre used in this work is 'Super 6T2' 
in the 'Laboratoire Leon Brilloin' in Saclay (France), where a resolution of 0.15°vertically 
and 0.1 "horizontally can be achieved. We have measured the intensity distribution of the 
as-grown crystals around the four antiferromagnetic Bragg reflections (|, =r, |), (\,\,\), 
(§' I' tO an d (t 1 ' h §)■ T ne P ea k splitting only occurs along one of the three symmetry 
allowed directions as shown for the (|, |, |) reflection in Fig. 1-a. Therefore, the modulated 
structure has a unique propagation vector k\ = (50 — 5) with 5 = 0.0064(1) corresponding 
to a period of 64nm. The elongated shape of the measured sattelites is due to the better 
resolution in the horizontal direction (along (10-1)) but also possibly because of a slight 
warping of the sample induced by the silver epoxy electrodes apposed on both sides of the 
crystals (corresponding to the (010) plane) for electric poling. 

The spin rotation plane can also be determined because the magnetic scattering amplitude 
depends on the component of magnetic moments perpendicular to the scattering vector. A 
quantitative analysis of the integrated intensities of 10 theta/two-theta magnetic reflections 
(see table I) allows us to conclude unambiguously that the moments lie in the plane defined 

— * — * 

by k\ = (50 — 5) and the polarisation vector P//[lll] (fig. 1-b). 

The structure refinement also confirms that the periodic structure is indeed a circular 
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cycloid with antiferromagnetic moments //p e = 4.11(15)/zb. Using a SDW model, or in- 
troducing a 20% ellipticity, deteriorates significantly the agreement factor of the fit. A 

— * 

consequence of the single k vector of the cycloid is that the crystal symmetry is lowered. 
Indeed, the ternary axis is lost and the average symmetry becomes monoclinic with the 
principal direction along /c=(110) [13]. 

No electric field effect on the magnetic order has ever been reported in bulk BiFeO^. Here, 
we analyse the effect of poling in the (010) direction perpendicular to the platelet. Because 
polarisation changes force charges to re-organise, the polarisation state of the sample can be 
monitored by measuring the current in the circuit. As the first coercive field was reached, 
we measured a significant decrease in the (|, |) neutron reflection intensity, indicating 
a redistribution of the average rhombohedral distortion. After reaching a multidomain 
state with < P >!=» 0, several neutron diffraction scans were performed. When trying to 
map precisely the intensity distribution around the antiferomagnetic Bragg positions, we 
found that the vertical resolution used for the crystal in its virgin state (in fig. 1) was 
not sufficient. Indeed, in the multi-domain state, ferroelastic distortions twin the crystal 
and complicate the diffraction patterns. In order to obtain a meaningful measurement, we 
had to reach 0.1 "of resolution in both horizontal and vertical directions, which pushes the 
experimental conditions to the limit of what can be done with these instruments. Figure 
2-b shows the (3D) reciprocal space mapping of the crystal. Yellow (111) type reflections 
are purely nuclear in origin while the red (|, |, |) are purely magnetic. (Ill) and (1-11) 
reflections are split along the long diagonals (dashed lines), which indicates the presence 
of two domains with different reticular distances. These are two rhombohedral twins with 
polarization axes along (111) and (1-11), roughly 50% — 50% in volume. The other (-111) 
and (11-1) reflections are also split, but along the (101) direction. This is due to a buckling 
of the crystal schematically shown in fig. 2-a, which slightly changes the angles fulfilling the 
Bragg conditions. This is fully consistent with polarised optical microscope images taken on 
similar crystals (fig. 2-a) indicating that the multi-domain state consists of stripe regions 
with two different polarisation directions. 

The purely antiferromagnetic peaks have been analysed in more details. The strongest 
(§' ~T> D reflection is shown in the zoomed region of fig. 2-b to be composed of four spots. 
These result from two simultaneous splits, one due to the ferroelectric distortion (already 
evidenced in the nuclear peaks) and one of magnetic origin. A projection of the zoomed area 
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is represented in figure 3 on with the expected reflections from P m and Pi-u domains are 
shown as green spots. The magnetic satellites are also indicated as black spots for the cycloid 
in the original (-101) direction and white spots for the other two symmetry allowed ones. The 
domains with the original Pm direction of polarisation lead to the pattern in the lower half 
of the figure, while those where the polarisation rotated by 71°(Pi_n) are in the upper half. 
In the latter domain, the expected satellites are not in a regular rhombohedral symmetry 
because they do not belong to the (111) diffraction plane of the figure. However, were they 
present, these satellites would still appear because the figure is a projection. Figure. 3 shows 
that in both domains, the splitting is only in the horizontal direction. Hence, the original 
(50 — 5) propagation direction of the cycloid in the virgin state was retained everywhere. 

The rotation planes of the AF vectors in the two domains can again be determined using 
the integrated intensities of the magnetic reflections (table II). These can be well accounted 
for by considering that 55% of the crystal volume has switched its polarisation by 71°, and 
brought with it the rotation plane of the Fe moments. Thus, in each domain, AF moments 

— * — * 

are rotating in the plane defined by k\ and P as represented in fig. 4. Hence, the electric 
field induced change of polarisation direction produces a spin flop of the antiferromagnetic 
sublattice. 

These measurements unambiguously demonstrate that the magnetic Pe 3+ moments are 
intimately linked to the polarisation vector. This negates the common belief that in bulk 
BiFe0 3 magneto-electric coupling must be weak because the cycloid cancels linear ME ef- 
fects [14, 15, 16, 17]. Although (M) = imposes a zero global linear ME effect in the bulk, 

— * — * 

the atomic coupling between M and P still exists. The underlying relevant microscopic 
mechanism is the (generalised) Dzyaloshinskii-Moriya (DM) [20] interaction which has re- 
cently been re-addressed starting from electronic Hamiltonians including spin-orbit coupling 
[18, 19]. Katsura et al. [18] describe in terms of spin currents the polarisation induced by a 

— * — * — * — * 

cycloidal spin arrangement, which can be written as P oc elj x (Si x Sj), with S^j the local 
spins and elj the unit vector connecting the two sites. The interaction of this polarisation 
with a coexisting internal polarisation produces a magneto-electric term in the total energy 
[18]: Edm = (P x elj). (Si x Sj). This ME interaction, which can also be obtained from 
symmetry considerations [12, 17], was held responsible for the cycloidal spin arrangement in 
BiFe0 3 [17]. This coupling energy induces the canting of Fe moments which exactly com- 
pensates for the loss in exchange energy (neglecting the anisotropy energy): E = —Ak 2 . A 
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ME energy density of — 3. 10 7 J/m 3 can be inferred from the value of the period of the cycloid 
and the exchange constant (A = 3. 10~ 6 J/m). Importantly, the coupling energy is zero when 

— * 

P is perpendicular to the local moments and maximum when it lies in the cycloid rotation 
plane. This explains the antiferromagnetic flop we observe when P changes direction. This 
also explains why the two crystals we measured had their cycloids in the same direction k\. 
Indeed, this minimises the components of the magnetic spins parallel to the depolarisation 
field (normal to the platelets surface), which lowers the cost in DM energy. 

When in thin film form, BiFe0 3 is a very different system because epitaxial strain sup- 
presses the cycloid and induces a weak magnetic moment [21]. Locally, the magnetic struc- 
ture consists of canted spins with angles changing sign between neighbours, which makes 
the moments add. If this magnetic configuration were to generate a local polarisation, its 
direction would alternate from site to site. Therefore, in order for a global polarisation to 
coexist with weak ferromagnetism, it is better if the spins are in a plane perpendicular to the 
polarisation direction, a configuration for which the DM based interactions are zero. This is 
exactly what is observed in BiFe0 3 films [5]. The magnitude of this ME coupling is more 
difficult to estimate than that in the bulk, but because it originates from the frustration of 
the DM interactions, it is likely to be weaker. Interestingly, canting angles are only about 
0.2°(the weak ferromagnetic moment being 0.02/is/atom [21]), to be compared with 2.25°in 
the bulk. This underlines the ME origin of the cycloid and also hints at a stronger coupling 

— * — * 

in the bulk since the interaction between P and M is directly linked to canting. 

In order to make a useful device with BiFe0 3 , we suggest to use the exchange bias inter- 
action between a thin ferromagnetic layer and a BiFe0 3 substrate, i.e. with its cycloid. It 
should be possible to vary electrically the exchange bias interaction using the antiferromag- 
netic flop observed here. Indeed, it is known that in conventional exchange bias systems, a 
non-compensated antiferromagnetic surface is not a prerequisite to obtain a large exchange 
field. Hence, it is likely that the cycloid may not significantly affect the bias, while optimising 
the coupling between the antiferromagnetic and ferroelectric orders. 
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and Daniil Khomskii as well as funding from the 'Agence Nationale de la Recherche' through 
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FIG. 1: (a) Neutron intensity around the (5, 2) Bragg reflection in the single domain state. 
The two diffraction satellites indicate that the cycloid is along the (10-1) direction, (b) Schematics 
of the magnetic configuration of antiferromagnetic vectors in the 64nm periodic circular cycloid. 
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FIG. 2: Mapping of the neutron intensity in reciprocal space. Two sets of splitting appear for the 
nuclear intensity (yellow spots) consistent with the presence of two ferroelastic domains (see (a)): 
one because of the presence of two rhombohedral distortions along [111] and [1-11], and the second 
because of a physical buckling of the crystal induced by the twinning. Magnetic peaks are further 
split because of the cycloids. Note that because the splitting is small, the scale has been magnified 
by a factor of 10 on each peak position. 
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FIG. 3: Neutron intensity around the (i, g) Bragg position in the multidomain state. Theo- 
retical positions are indicated by the black and white spots. Diffraction satellites are visible in the 
0°(bottom half) and 71°(top half) domains of polarisation. The difference in vertical spot shape 
likely originates from the position of reversed domains at opposite ends of the sample because of 
a prefered nucleation near the the edges. Any warping of the sample splits the new (1-11) peak 
while recovering an improved resolution for the original domain located near the center. 
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FIG. 4: Schematics of the planes of spin rotations and cycloids ki vector for the two polarisation 
domains separated by a domain wall (in grey). 
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TABLE I: Intensity measured around the magnetic Bragg positions compared to that expected for 
a cycloid with magnetic vectors in the different allowed (P, k) planes. 
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TABLE II: Intensity measured around the magnetic Bragg positions compared to that expected for 
a cycloid with magnetic vectors in the different allowed planes. Calculated values are obtained with 
55% of domains having switched their polarisation by 71°and kept the same propagation vector. 
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